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We discuss the description of quantum magnetization in the super paramagnetic compound Fes 
using a generalization of the Lipkin-Meshkov-Glick Hamiltonian. We study the variation of the 
energy spectra and of the wave-functions as functions of the intensity of an external magnetic field 
along the three magnetic anisotropy axes. 
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I. INTRODUCTION 

Single Molecule Magnets (SMM) are polynuclear co- 
ordination compounds that show slow relaxation of the 
magnetization [5j . Among their features we can find high 
electronic spin ground state, a ferromagnetic or anti- 
ferromagnetic coupling due to super exchange interac- 
tions, i. e. mediated by diamagnetic bridge ligands; and 
an organic ligands shield protecting their magnetization 
from the environment Q. Beyond the unusually high 
spin value, the main feature of the SMMs is the ten- 
dency of each molecule to remain magnetized after ax 
plying and then removing an external magnetic field 
The relaxation rate is described by the equation: 



t = tq exp 



U 



(i) 



where r is the relaxation time, To is the relaxation time at 
infinite temperature, U is the magnetic anisotropy con- 
stant, ks is the Boltzmann constant and T is tempera- 
ture. The cluster Mni2 Ac has exhibited a relaxation time 
of about two months at liquid helium temperature 
The shown qualities can be explained with the giant spin 
model. According to it, the interactions between para- 
magnetic ions within the molecular assemble produce an 
effective magnetization of the whole molecule [2]. Be- 
cause of this phenomenon, this kind of systems can be 
described by the Zero Field Splitting Hamiltonian: 
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where S is the spin operator, D is known as the 
anisotropy tensor, D and E are the axial and rhom- 
bic anisotropy parameters respectively, and the Si(i = 
x,y, z) are SU(2) algebra operators acting over the spin 
states at each Cartesian direction 12]. Such Hamiltonian 
refers to the magnetic anisotropy which is related to the 
molecular geometry, making the two terms included in 
the ZFS Hamiltonian dependent on the symmetry of the 
molecule |2j . The interaction of the magnetization with 
an external magnetic field gives rise to the Zeeman Effect, 
which is given by a Hamiltonian of the kind: 



(3) 



where g is the electron gyromagnetic constant, is the 
Bohr magneton and Bq is the external magnetic field 



12]. Amongst all the contributions to the general de- 
scription of the electronic spin, the two terms of the spin 
Hamiltonian mentioned above lie in the energy range of 
to 1 cm _1 )[l2j], making them the most relevant for the 
magnetization analysis in SMM. The sum of both terms 
is written as: 

H = D (Sl - |S 2 ) + \E (Si + S 2 _ ) + gfi B S • B (4) 



This formulation corresponds to a direct generalization 
of the model developed by H. J. Lipkin, A. J. Glick and 
N. Meshkov in 1965 to study the many body problem in 
nuclear physics (lfjj . This model and its mathematical 
and physical implications had been widely studied aside 
from its similarity with the molecular expression [9[ . The 
rising of chemical species like the SMMs produced the 
recognition of the relation between both formulations. 

As a consequence of their particularities, the SMM 
are quite good candidates for a huge variety of appli- 
cations, including the image magnetic resonance as con- 
trast agents @ , as information storing devices and as pro- 
cessors in quantum computing (l3j | . An implementation 
of the Grover's search algorithm with pulse sequences 
in Electron Paramagnetic Resonance experiments using 
SMM has been proposed [III- Entanglement of SMM 
with cavities is a candidate for W states generation [lij ]. 
Molecular dimmers also exhibit quantum entanglement 
Q. The observation of quantum interference associated 
with tunnelling trajectories between states of different 
total spin length in a dimeric molecular nanomagnet pro- 
vides clear evidence for quantum-mechanical superpo- 
sitions involving entangled states shared between both 
halves of wheel-shaped molecule 15[ . 



Fe 8 



The octanuclear SMM [Fe 8 2 (OH)i2(tacn) 6 ] 8+ , where 
tacn is the organic ligand 1,4,7-triazacyclenonane, known 
simply as Fes (figure 1) is of particular interest due to 
its approximate D2 symmetry, allowing all the terms in 
expression Q to be active in its description [3(. 

The metallic cluster is formed by eight Fe 3+ ions in 
high spin configuration (s = 5/2) bridged by oxo and 
hydroxo ligands. The ground state with S = 10 is de- 
scribed considering that six of the eight ions have their 
spins aligned ferromagnetically, while the remaining two, 
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nian ((4]). Situations in which the external magnetic field 
is aligned with different directions referred to the axes es- 
tablished by the molecular magnetization are discussed. 
A slightly improved description of the magnetization of 
this particular molecule can be obtained by the addition 
of Stevens operators involving fourth order spin operators 
Q . There are also experimental evidences of a small but 
measurable contribution of the S = 9 state when working 




Figure 1: Representation of the Feg molecule [131]. Color code: 
Fe, orange; O, red; H, white; N, purple and CH2, gray. 



Table I: Values of the anisotropy constants reported in liter- 
ature. 

Author Experimental Technique D [K] E [K] 



Barra [3] 
Caciuffo [1] 
Hill [7] 



EPR 
INS 
EPR 



-0.275 -0.046 
-0.292 -0.047 
-0.292 -0.046 



/hfill 

at temperatures higher than those studied in this work 
'16]. These two aspects are not included in the present 
contribution. 



the ones surrounded by the greater number of Fe 3+ 
ions at second neighboring; are aligned antiferromagnet- 
ically in the arrange |3| . The anisotropy parameters had 
been experimentally determinated using techniques such 
as electron paramagnetic resonance (EPR) and inelas- 
tic neutron scattering (INS) They are listed in Table U 
where it can be appreciated that the value of the axial 
anisotropy coefficient is negative, which would imply that 
the molecule ground state could have to the largest value 
of the spin projection (Ms) 0- Nonetheless, the presence 
of the rhombic anisotropy coefficient produces eigenfunc- 
tions of the Hamiltonian that can be expressed as linear 
combinations of the functions given by the spin projec- 
tions. The quadratic characteristic of the operators make 
that the superpositions involve only same parity projec- 
tions m. 

The present work gives a detailed study of the exact 
solution of the Schrddinger equation with the Hamilto- 



II. RESULTS AND DISCUSSION 

The results presented here were obtained through 
numerical calculations employing MATLAB ©. The 
anisotropy constants used were taken from Hill et. al. 
(Table H 

A. Wavefunctions 

The energy levels and their associated wave-functions 
are the eigenvalues and eigenvectors of Hamiltonian Q : 

H fa) = £i |V><) , (5) 
which matrix elements are: 



(Ms\H\M' s ) = D [M 2 S -\S(S+ 1 )] + gi i B B Q cos 9M S 

+ ytiBBg smOe-V y/S (5 + 1) - M s (M s + l )S Ms , M > s+1 
+ \gn B B smBe^y/S (S + 1) - M s (M s - l)5 Ms ,M's-i 

+ \eJ [S (S + l)f - [2[S (S + l)f - M 2 S - 2M S ) (Ms + l) 2 <5 Ms , M ' s+ 2 
+ \E\J [S (S + l)] 2 - [2[S (S + l)] 2 - Ml + 2Ms\ (M s - l) 2 6 Ms M's-2, (6) 

I 



where <j) is the angle in the xy plane, 9 is the angle with magnetic field. The eigenfunctions in ([5]) are of the form: 
the z axis and Bq is the total intensity of the external 

s 

\i>i)= J] Ci >Ms \M s ) (7) 



M s =-S 
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Figure 2: Energy spectrum of the Fes spin states varying the intensity of the magnetic field aligned with the easy molecular 
axis. Colors identify the main Ms component at zero field. Right: Calculation assuming biaxial anisotropy. Left: Calculation 
assuming only uniaxial anisotropy. 



When the field B is aligned with the z axis (easy axis) 
all the terms in the linear combination have the same 
parity, i. e. they correspond to only even or odd values 
of Ms- Such situation generates a spectrum as the one 
shown in figure 2-right. 

With a merely illustrative aim, we show also the 
hypothetical spectrum in which there is not rhombic 
anisotropy, i. e. E = (figure 2-left). In this case it is 
observed that the dependence of each energy level with 
the magnetic field is perfectly linear, in agreement with 
what is commonly observed for the Zeeman Effect, while 
for the real one the dependence is polynomial, conse- 
quence of the variation, with the magnetic field, of the 
coefficients multiplying the spin projections on each lin- 
ear combination.. The former produces interesting situ- 
ations not observed with pure axial anisotropy. 

In the realistic situation at high field an approximately 
linear behavior is observed. On the other hand, with field 
values smaller than 1 T, when the energy lines between 
two levels get close, not always happens a crossing be- 
tween them. This situation originates the coexistence of 
real crossing (two levels share the same energy at cer- 
tain value of the magnetic field) and avoided crossing 
(the curves of two levels get closer up to a point where 
the slope sign changes, avoiding them to coincide at a 
specific coordinate) (figure 3). In a real crossing the co- 
efficients multiplying the spin projection on each level 
remain practically unchanged in the vicinity of the cross- 
ing. On the contrary in an avoided crossing the values of 
the coefficients are exchanged between the levels taking 
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Figure 3: Avoided crossing between the 11 th and the 12' 
excited spin states (steady lines). The asymptotic behavior 
reveals a correspondence between the real crossing with spin 
projections -3 and 7 (dashed lines). 



part in the rapprochement. For the Hamiltonian (TjTJ with 
the external field aligned with the easy axis the avoided 
crossings can take place only between levels having the 
same parity of the main spin projection. 

It is appreciable that the higher energy levels have a 
behavior very far from linear at low field, producing very 
pronounced avoided crossings. It is a signature that, in 
this region, the mixing between spin projections is large 
in the states associated with these energy levels. 
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Figure 4: Probabilities of finding the different spin projections along the easy axis obtained from the wave function of the 4 th 
excited state at different ranges of magnetic field intensity. The average expectation value of S z for each magnetic field region 
are given in the insets, the consistency between this value and the dominant components on each case is noticeable. 



In figure 4 the average spin projections over the easy 
axis are included in the insets, and compared with the 
probabilities of finding this particular spin projection in 
the wave function of the 4 th excited state, at different 
magnetic field ranges. Those ranges were splitted follow- 
ing the notorious predominance of one of the spin projec- 
tions over any other, predominance that changes at the 
values of magnetic field where crossings take place. This 
value is in good match with the predominant component. 

In the cases in which there are components of the field 
along the other directions, perpendicular to the easy axis, 
the linear operators S x and S y start to be relevant and 
the components of the wave function have no parity re- 
strictions. This generates significant modifications on the 
energy spectrum, having situations in which the tendency 
of mixture is changed and even their relation with the 
energy at low field is inverted (figures 5. a and 5.b re- 
spectively). Moreover, a maximum mixture is obtained 
when the external magnetic field is aligned with non of 
the internal axes (figure 5.c). 

It can be noticed that with some alignments, the cross- 
ing zone keeps regions with some sort of reticular struc- 
ture where the avoided crossings coexist with the real 
ones, the persistency of such substructures is due to the 
low relative value of E regarding to D. 

As an example of the importance of the effect of the 
level mixture, figure 6 shows the variation in the squared 
coefficients of the spin projections taking part in the wave 
function of the 4 th excited state as a function of the ex- 
ternal magnetic field aligned with the hard axis. It can 



be observed that, at variance from the alignment with 
the easy axis, there are no field intensity ranges in which 
the predominance of a single component over any other 
is recognizable. An interesting fact is that projections 
with the same absolute value produce the same curve. 
All these contributions produce that the expected value 
(S z ) remains around zero over all the range of magnetic 
field. 



III. CONCLUSIONS 

We have analyzed the wavefunctions describing the 
magnetization in the molecular cluster Fes- We have 
shown that, with a given molecular symmetry, the be- 
havior of the wavefunction and of the energy levels in 
the presence of an external magnetic field holds a strong 
relationship with the alignment between that field and 
the molecular axes. In this case the spin operators act 
like actual geometric components of the quantum magne- 
tization within the molecule, so it is easy to assign an op- 
erator with a perfectly defined axis, producing the coor- 
dinate system to rely on the molecule's orientation. The 
molecular magnetic anisotropy has, as a consequence, the 
mixing of the states in the spin basis. This mixing also 
depends on the alignment with an external non-zero field. 
This prominent feature which can be exploded for quan- 
tum algorithms. Phase transitions and entaglement in 
collective spin systems have been studied in a theoret- 
ical framework previously (see Ref. Q and references 
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Figure 5: Energy spectra of Feg as function of the external magnetic field. Colors identify the corresponding values of the 
dominant Ms component at zero field, a) Field aligned with the hard (x) axis, b) Field aligned with the intermediate (y) axis, 
c) Field at 45° of the three axes. 
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Figure 6: Variation of the square of the coefficients of the spin 



projections contained in the wave function of the 4 
state as a function of the magnetic field. 



excited 



therein). In the next stages of the research we will inves- 
tigate specfically applications of SMM to quantum infor- 
mation, in particular entanglement and quantum phase 
transitions. 
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